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ON EXTRACTION OF THE TOTAL PHOTOABSORPTION CROSS SECTION ON THE 
NEUTRON FROM DATA ON THE DEUTERON 
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An improved procedure is suggested for finding the total photoabsorption cross section on the 
neutron from data on the deuteron at energies < 1.5 GeV. It includes unfolding of smearing 
effects caused by Fermi motion of nucleons in the deuteron and also takes into account non- 
additive contributions to the deuteron cross section due to final-state interactions of particles in 
single and double pion photoproduction. This procedure is applied to analysis of existing data. 



Introduction 
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This work was motivated by recent preliminary results from the GRAAL experiment on the total 
photoabsorption cross section off protons and deuterons at photon energies ft) = 700—1500 MeV 
Q3-31 and their implications for the neutron. An intriguing feature of the new data is that they 
indicate an approximately equal and big strength of photoexcitation of the nucleon i<i5(1680) 
\ resonance off both the proton and neutron (as seen, in particular, in Fig. 5 in Ref. flU). Meanwhile 

this strength was found small for the neutron in many previous studies (see, e.g., (HO). Particle 
Data Group [7] quotes the following branching ratios of /V* = Fis(1680) to yN: 

*t '• 

73: Br(N*^yp) = 0.21-0.32%, 



Br(/V*->yn) = 0.021-0.046%. (1) 



Irrespectively on whether the old or new data are correct, it seems timely to (re)consider procedure 
commonly used to find cross sections off the neutron from the deuteron data. 

This procedure was described in detail by the Daresbury group [6] who performed measure- 
ments of the total photoabsorption cross sections o p [5| and a c i [6| at energies between 0.265 and 
4.215 GeV. In the nucleon resonance energy region they made an Ansats that 



o d (co)=F(o))[o p ((o) + G n (a})]. (2) 

Here the factor of F(co) was introduced in order to take into account smearing effects due to Fermi 
motion of nucleons in the deuteron. This factor was found by numerical integration of the proton 
cross sections using known momentum distribution of nucleons in the deuteron and then equally 
applied to the neutron. Finally, the neutron cross section was found, point by point, with the step 
of 25 MeV, from the corresponding deuteron cross section at the same energy using Eq. ©. 

An evident drawback of the Ansatz © is that smearing effects are assumed to be the same for 
the proton and neutron, what cannot be true in case the energy dependencies of o p (co) and o n ((Q) 
are different. 

The second problem is that smearing of the cross section makes it impossible to relate individ- 
ual nucleon cross sections C)v(ft)) with o^(ft>) at the same energy and thus to apply the point-by- 
point procedure. Instead, some average of O)v(ft)) over a finite energy interval can only be found. 
In other words, a justified unsmearing procedure should be applied there. 

The third point is that non-additive corrections related mostly with final state interactions have 
been neglected in Eq. ©. Brodsky and Pumplin [8] estimated these corrections at high energies 
(ft) > 2 GeV) assuming that high-energy photoproduction on the nucleon is dominated by diffrac- 
tive photoproduction of vector mesons (p, ft), 0) which then interact with the second nucleon. 
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Such corrections have been included in the analysis of high-energy part of the Daresbury data [|6] 
(as well as in studies of photoabsorption off protons and deuterons at energies 20-40 GeV [9|). 
At lower energies, including energies of GRAAL, the corrections related with vector meson pro- 
duction are small. Nevertheless, other photoproduction channels still might be important. This is 
indeed the case as explained below. To our knowledge, no estimates of the non-additive corrections 
to Eq. ([21) have been yet done at energies of the GRAAL experiment. 

In this work we improve the procedure of [16] in all the above three lines. 

Fermi smearing (folding) 

We begin with rewriting Eq. © more accurately as 

o d ((o)=F[o p ((o) + o n ((o)}+AOp n ((o). (3) 

Here F is a linear integral operator that smears individual nucleon cross sections in accordance 
with Fermi motion of nucleons in the deuteron; Ao pn is a non-additive correction to be discussed 
later. The first two terms in Eq. © arise from diagrams of impulse approximation (like those in 
Fig. [T]) when interference effects are omitted. We neglect here off-shell effects for intermediate 
nucleons N because the binding energy of nucleons in the deuteron is rather small (2.2 MeV). 
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Figure 1 : Diagrams of impulse approximation for yd — > nNN. Antisymmetrization over N\ and 
N2 is not shown. 

A simple analysis of diagrams of impulse approximation shows IfTOl that the smearing operator, 
in nonrelativistic approximation over nucleons in the deuteron, is reduced to 

W{p z ) — o N {af E )dp z . (4) 



of s =( (l-E±) (5) 



Here 

m eff = m( 1 - 

M, 

is the effective (Doppler shifted) energy for the moving intermediate nucleon N of the mass M pro- 
vided its longitudinal (along the photon beam) momentum is equal to p z . W(p z ) is the longitudinal 
momentum distribution of nucleons in the deuteron, 



W( Pz )= I \¥(P)\ 2 ^, (6) 



d 2 p 

and the factor co eff / £0 takes into account a change in the photon flux seen by the moving nucleon. 
As in Ref. [61, we use in the following a simplified deuteron wave function (Hulthen limo . 



W ( r ) = -(e- ar -e- br ), I | y{r)^A%r L dr = 1, (7) 



r Jo 
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with a = 45.7 MeV/c, b = 260 MeV/c and k 2 = ab(a + b)/[2n(a-b) 2 ] = 12.588 MeV/c. In the 
p-space 

so that the function W(p z ) is 

<^> = ^(i + 5-^?f)> /»W*=i, CO 

where A = a 2 + p 2 and B = b 2 + p 2 . This function is shown in Fig. [2] together with a distribution 
obtained with a realistic (CD-Bonn) wave function [12]. In actual calculations we cut off mo- 
menta \p z \ > p cut = 200 MeV/c where W(p z ) becomes quite small and the momentum p z remains 
nonrelativistic. 




-0.2 -0.1 0.1 0.2 

P z (GeV) 

Figure 2: Distribution of the longitudinal momentum in the deuteron. Solid and dashed lines: 
Hulthen and CD-Bonn wave functions. 



The Hulthen distribution for W(p z ) gives the following average longitudinal momentum of 
nucleons in the deuteron: 

(pl) l / 2 = 53.9 MeV/c (10) 

(it is 54.9 MeV/c for the CD-Bonn wave function). It also gives the following spread for the 
effective photon energy seen by the moving nucleon: 

/r,2\l/2 

Aco ee = co^j— = 0.057w. (11) 

In other words, this value characterizes the "energy resolution of the deuteron" as a "spectral 
measuring device" for the neutron. For CO ~ 1 GeV only an average of the nucleon cross section 
over the range ~ ±60 MeV can be inferred from the deuteron data. Determination of O n ((o) with 
the step of 25 MeV done in |[6l cannot be physically justified. 

Unfolding 

It is well known that the unfolding problem, i.e. solving the Fredholm integral equation © for the 
unknown "unsmeared deuteron cross section" a(co) = o p (co) + o n ((Q), cannot be solved without 
further assumptions on properties of the solution o(co). In particular, it is not possible to restore 
fast fluctuations in o(co) at the energy scale < Act) . To proceed, we make therefore a physically 
sound assumption that both the cross sections (7 p (oj) and <7„(co) can be approximated with a sum 
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of a few Breit-Wigner resonances (having fixed known standard masses and widths but unknown 
amplitudes, probably different for p and n) plus a smooth background. Thus we write 

<J(«)=£X ( - fi(co) (12) 

where fi(co) is the basis of the expansion, i.e. either Breit-Wigner distributions or smooth functions 
of the total energy ^fs. We borrow specific forms of the functions fi{(0) from Ref. [6J, Eqs. (11) 
and below. Then unknown coefficients Xj are determined from the fit of Fo(co) to experimental 
data on o c t(oo) (at this point we assume that the correction Ao pn is already calculated). 

A knowledge of X,-, with errorbars 8Xi determined in the fit, can be directly converted to the 
knowledge of 0(00), also with errorbars. In particular, writing fluctuations in the determined value 
of o(co) as 

80(00) =^8Xifi((0), (13) 

i 

we have 

8a 2 ((o) = Y,SX i 8X j f i (oo)f j (o)) (14) 

ij 

and 

(5(7 2 (£0))=£C IJ / i (£0)/ ; («), (15) 

ij 

where 

C ij = (8X l 8X J ) (16) 

is a standard covariance matrix of errors determined in the fit of Xi. 

In this way the extracted unfolded cross section a (ft)) can be shown as a smooth curve (corre- 
sponding to the central values of X,) surrounded with a band of the half-width given by Eq. <TT3T > 
which represents errors in the cross section. 



Nonadditive corrections 

The term Ao pn ((o) in Eq. ([3]) takes into account various effects violating additivity of the photoab- 
sorption cross sections on individual nucleons. Among them: 

- interference of diagrams of photoproduction off proton and neutron, Fig. \T\ leading to iden- 
tical final states; the Fermi statistics of the emitted nucleons (antisymmetrization) leading to the 
so-called Pauli blocking, 

- interaction between emitted particles (final state interaction, FSI) including both interaction 
of unbound nucleons and binding of nucleons (formation of the deuteron in the final state), in- 
teraction of pions (or other particles), produced on one nucleon, with the second nucleon in the 
deuteron, 

- absorption of pions (and the presence of processes such as the deuteron photodisintegration, 
without pions in the final state). 

Now we briefly discuss all these effects starting with the reaction of single-pion photoproduc- 
tion, yd — > nNN, considered in the model that includes diagrams of impulse approximation (Fig.[T) 
and the final state AW and TtN interaction to one loop (Fig. [3). Formalism and the main building 
blocks of this model that was previously used in the energy region of the A(1232) resonance can 
be found elsewhere [|T3llT4Tl . Generally, the model works well for the channel yd — > ll~pp in the 
A(1232) region but not so well for yd — > n°pn, see Fig. 0] Reasons for the discrepancy are not 
clear but other authors get similar results and also cannot describe the data (see, e.g., [[T7T0 . We 
will not use the model for energies too close to the A(1232) region. 
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Figure 3: Diagrams with the final state AW and 7tN interaction (to one loop) for yd — > TtNN. 
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Figure 4: Model H13U141 predictions for yd — > % pp (left) and yd — > n°pn (right) in the region of 
the A( 1232). 



In the present calculation that covers higher energies, "elementary" amplitudes of yN —> TtN 
are taken from the MAID analysis [fT5l (with a proper off-shell extrapolation); those for NN — > NN 
are taken from the analysis of SAID [16] (again with an off-shell extrapolation). In the following 
plots we show obtained results for Ao pn (co) in different isotopic channels. 



1. Interference contributions from diagrams of impulse approximation for 7<i — y irNN 
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Figure 5: Contribution to Ao pn due to interference of diagrams in Fig. [T]of impulse approximation 
for yd nNN. 
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2. NN FSI interaction in jd — > ttNN and *yd — > nd 

We put here AW FSI contributions for the continuous and bound states together because there is a 
tendency for their cancelation that can be traced to the unitarity (closure). The matter is that the 
AW interaction in the continuous spectrum can be thought as a replacement of the plane AW wave 
in the reaction amplitude of the plane- wave impulse approximation, 

T pwlA (E NN ) = (NN\T(YN-+nN)\d), (17) 

with the distorted AW wave in the reaction amplitude of the distorted- wave impulse approximation, 

T DW1A (E NN ) = (^-\NN)\T(rN^nN)\d). (18) 

Here we explicitly indicate the energy of the NN state. Also, the coherent amplitude, with the final 
bound AW system, is 

T coh (E d ) = (d\T(yN^nN)\d). (19) 

Owing to the closure, i.e. a completeness of eigen states of the free AW Hamiltonian as well as 
those of a Hamiltonian with AW interaction, 

1= £ |AW)(AW|= £ \y(-\NN))(y(-\NN)\ + J»R (20) 

NN,E NN NN,E NN d 

the square of the PWIA off-shell amplitude integrated over all possible AW states, irrespectively to 
their energies, exactly coincides with the square of the DWIA off-shell amplitude (also integrated 
over all possible states) plus the square of the coherent amplitude. In case when a subset of AW 
states of certain energies is only considered, as in the case of finding cross sections at a certain en- 
ergy, the coincidence of |r PWIA | 2 with |j DWIA | 2 _|_ |ycoh|2 - ys not strictly valid, however a tendency 
to have a compensation between the coherent contribution to the cross section and a decrease in 
the DWIA cross section still remains. 

An illustration of this general tendency can be found in Fig. [6] where the negative AW-FSI 
contribution to yd — > n°pn is close in the magnitude to the positive coherent contribution to yd — > 
n°d (see dotted curves). 




iC - E i i i i i i i ' ■ h t — ■ — ■ — ■ — ■ — 1 — ■ " ' i -1 
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Figure 6: Left: Contribution to Ao pn due to final state NN interaction in yd — > nNN. Right: 
Contribution to Ao pn from yd — > n°d and yd — >■ %%d. 
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3. NN FSI interaction in jd — > ttttNN and jd — > irnd 

Consideration of the reactions yd — > tzkNN and yd — > KKd is similar but more involved owing 
to a more complicated structure of the elementary yN — > TtnN amplitude. We rely here on results 
obtained by Fix and Arenhovel HI 9 11201 from which we infer contributions to Aa pn shown in Figs. [6] 
(the right panel) and [7] Again we see an essential partial cancelation between yd — > n + iz~d and 
AW -FSI effects in yd — > lz + K~pn. 

5 
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Figure 7: Final state AW interaction in yd — > nnNN. 




4. Other small contributions and the net result for Acr pn 

We do not show contributions to Ao pn from tzN FSI in yd — > nNN (found in the described model) 
and contributions from the deuteron photodisintegration, yd — > pn (it can be directly found from 
experimental data of CLAS [18J) because they are rather small with the except for energies close 
to the A(1232) resonance region. We can anticipate that Ao pn is not affected by rj meson pho- 
toproduction because r]N interaction is weaker than that of TlW and because effects of AW FSI 
interaction in the continuum and in the bound state are again nearly canceled. 

Taking all contributions together, we arrive at the total value of Ao pn shown in Fig. [8] which is 
the main result of this section. In spite of quite a few pieces of order 10 /lb, the sum of all contri- 
butions to Ao pn is found surprisingly small, so that our improvement to the unfolding procedure is 
mainly reduced to a refinement in solving the integral equation. 
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Figure 8: Total value of Ao pn . 
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Extraction of the photoabsorption cross section on the neutron 

Known now all ingredients of Eq. ©, we can fit experimental data, determine the unsmeared 
deuteron cross section <j p + o n and then find the neutron cross section a n . We illustrate this proce- 
dure using Daresbury data [|5J[6) for the proton and the deuteron. 

Figure |9] (the left panel) shows a smooth fit (the curve labeled "tot") with Eq. (fT2l) to the 
experimental proton data and the result of its smearing with the smearing operator F. Separately 
shown is the contribution of resonances (and its smearing) and a smooth background. At the 
right panel of Fig. [9] a fitting curve is shown that, after smearing and adding Ao pn , comes through 
experimental data points (the curve labeled "totF"). 




Figure 9: Daresbury data for the proton (left) [5] and the deuteron (right) fl6l, their fit and smearing. 

From this fit the neutron cross section can be found as a difference, see Fig.[lOl(the left panel). 
In a similar way the neutron cross section can be found from Mainz data ll2TTl . Our results are 
shown in Fig. [10] (the right panel). Bands indicate errors in the found neutron cross sections there. 



Conclusions 

An improved procedure of extracting the total photoabsorption cross section on the neutron from 
data on the deuteron is proposed. It involves a more correct treatment of folding/unfolding of the 
Fermi smearing of individual nucleon contributions. 

Non-additive corrections are evaluated at medium energies where VMD does not yet work. 
They are relatively small in total but they might be more important in analyses of partial channels 
of photoabsorption. 

We hope that the obtained results will be useful for interpretation of the GRAAL data and 
future experiments. 
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Figure 10: Extraction of the neutron cross section o n from the deuteron data (Daresbury [6 | and 
Mainz ETTO . Original values of o„ from the Daresbury experiment are shown with open circles. 
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